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We present a Linearized Coupled Cluster (LCC) correction based on an Antisymmetric Product 
of 1 reference orbital Geminals (APlroG) reference state. In our LCC ansatz, the cluster operator 
is restricted to double and to single and double excitations as in standard single-reference CC 
theory. The performance of the APlroG-LCC models is tested for the dissociation of diatomic 
molecules (C 2 and F 2 ), spectroscopic constants of the uranyl cation (UOj" 1 "), and the symmetric 
dissociation of the H 50 hydrogen chain. Our study indicates that an LCC correction based on 
an APlroG reference function is more robust and reliable than corrections based on perturbation 
theory, yielding spectroscopic constants that are in very good agreement with theoretical reference 
data. 


I. INTRODUCTION 

In the past decade, quantum-chemical studies prove 
helpful in understanding and predicting chemical phe¬ 
nomena of small and medium-sized molecules ID 121- 
Routine applications to larger molecular systems are, 
however, hampered as conventional quantum chemistry 
methods are either too expensive or too approximate to 
guarantee reliable results. These problems are partic¬ 
ularly severe for molecules where strong electron cor¬ 
relation becomes important. Examples for strongly- 
correlated systems are, for instance, radicals, transition 
metals, and actinide compounds. Furthermore, well- 
established approaches for strong correlation scale poorly, 
often factorially, with system size. These drawbacks mo¬ 
tivate the development of new, unconventional electron 
correlation methods for strongly-correlated many-body 
systems that represent computationally cheaper alterna¬ 
tives to standard methods. 

One active field of research for strongly-correlated 
systems focuses on the development of wavefunction 
methods that use two-electron functions (geminals) to 
model the correlated motion of electrons [Ml- The 
most popular geminal-based approaches are the Antisym¬ 
metric Product of Strongly orthogonal Geminals 00- 
m (APSG), the Antisymmetrized Geminal Power [T5b 
115] (which is a special case of projected Hartree-Fock- 
Bogoliubov |19]). the Antisymmetric Product of Inter¬ 
acting Geminals 0 EUH5T] (APIG), Generalized Valence 
Bond 0 I52H55] (GVB), and the Antisymmetric Prod¬ 
uct of 1-reference orbital Geminals (APlroG) [36il38| . 
Specifically, APlroG provides a very good approxima¬ 
tion to the doubly occupied configuration interaction [32 
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(DOCI) wavefunction, at mean-field computational cost. 
The APlroG geminal creation operator reads 


b] = ajat 


E 


c>la|, 


(1) 


where aj,, are the standard fermionic creation opera¬ 
tors for spin-up ( p ) and spin-down electrons ( p ), c“ are 
the geminal coefficients, and the summation runs over all 
virtual orbitals. Specifically, the APlroG geminal coeffi¬ 
cient matrix has the form 


CAPlroG = 


( 1 ... 

0 0 

r p+1 

C 1 

r p +2 . 

c l 

■ cf\ 

0 1 

... 0 

c p+1 

c 2 

r p + 2 . 
l 2 

• cf 

VO i 

... 1 

c p+1 

Cp 

„P+2 . 
Cp 

• 4 ) 


where K denotes the number of basis functions, P 
the number of electron pairs, and the left sub-block of 
CAPlroG entails some reference determinant. The gem¬ 
inal matrix connects each geminal with the underlying 
one-particle basis functions. We should note that if we 
impose specific restrictions on the structure of the above 
matrix, we can deduce different geminal models 0. 

The electronic wavefunction is written as a prod¬ 
uct of geminal creation operators for all electron pairs 
P acting on the vacuum state, | APlroG) = Ilf ^*1)- 
Unique among geminal methods, the APlroG wavefunc¬ 
tion ansatz can be rewritten in terms of one-particle func¬ 
tions as a fully general pair-Coupled-Cluster-Doubles [40] 
(pCCD) wavefunction, 


|APlroG) = exp ( ^ 



Ci a a a a™ 


a i a i |$0>, (3) 


where |<I>o) is some independent-particle wavefunction 
(for instance the Hartree-Fock (HF) determinant). The 
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exponential ansatz of APlroG (c/. eq. @) ensures size- 
extensivity of the model. However, to recover size- 
consistency, we have to optimize the one-particle ba¬ 
sis functions. This can be done in a fully variational 
manner [23135], analogous to orbital-optimized Coupled 
Cluster [41], or using non-variational orbital optimiza¬ 
tion techniques |42j, 43]. A number of numerical studies 
on systems with strongly correlated electrons showed the 
superiority of the variational orbital optimization proce¬ 
dure over the latter ones [43]. We should note that due 
to the four-index transformation of the electron repul¬ 
sion integrals, the computational scaling of the orbital- 
optimized APlroG model deteriorates to O (A" 5 ) [5T j . Al¬ 
though restricted orbital-optimized APlroG is limited to 
close-shell systems, it has already proven to be a reli¬ 
able method for modeling strong electron correlation ef¬ 
fects in quasi-degenerate systems [331421143], single and 
multiple bond-breaking processes [44] [45], and actinide 
chemistry [46] , 


As all other geminal models, APlroG misses a large 
fraction of weak (dynamical) electron correlation effects. 
To address this problem and account for weak elec¬ 
tron correlation effects in the geminal reference wave- 
function, various a posteriori corrections have been pro¬ 
posed. These include models based on single- and multi¬ 
reference Perturbation Theory E3 S3®, Extended 
Random Phase Approximation [511 f52j . (Linearized) 
Coupled Cluster theory [33 [35], and Density Functional 
Theory [34] [53- In the case of APlroG, dynamical corre¬ 
lation was included using Perturbation Theory [48] and 
Density Functional Theory [54, 53 . Recent studies on 
diatomic molecules, however, point out numerical insta¬ 
bilities and failures of the proposed Perturbation Theory 
corrections [45] . This motivates the development of dif¬ 
ferent, ideally more robust dynamical correlation models 
for APlroG. A reliable way to account for dynamical 
correlation effects a posteriori is to use a multi-reference 
Linearized Coupled Cluster (LCC) correction. Recently, 
Zoboki et al. presented an LCC correction based on an 
APSG reference function and demonstrated the good per¬ 
formance of the APSG-MRLCC approach. Their findings 
encouraged us to develop an LCC correction based on an 
APlroG reference state. 


This work is organized as follows. In section [IT] we will 
discuss two different LCC corrections for APlroG. Their 
performance is compared by studying some well-known 
problems in quantum chemistry that require a balanced 
treatment of dynamical and strong electron correlation 
effects: the dissociation of C 2 and F 2 , the symmetric 
dissociation of H 50 , and spectroscopic constants of the 
V0 2 + molecule. Computational details are presented in 
section ITT] while numerical results are discussed in sec¬ 
tion [IV Finally, we conclude in section [V] 


II. LCC THEORY WITH AN AP1ROG 
REFERENCE FUNCTION 

In this work, dynamical correlation effects are built in 
the electronic wavefunction a posteriori using an expo¬ 
nential Coupled Cluster ansatz, 

|'P) = exp(T)|APlroG), (4) 

where T = ^ v t v T v is a general cluster operator. The 
corresponding time-independent Schrodinger equation 
reads 

H exp (f) | APlroG) = A exp(T)| APlroG) (5) 

Multiplying from the left by exp(—T) and truncating 
the Baker-Campbell-Hausdorff expansion after the sec¬ 
ond term, 

exp(-f)Hexp(f)«H+[H,f], (6) 

we arrive at the Linearized Coupled Cluster Schrodinger 
equation 

(H+[H,T\) | APlroG) = E\ APlroG). (7) 

To obtain the cluster amplitudes i„, we multiply from 
left by (iz| 

{v\{H + [H,T]) | APlroG) = 0, (8) 

where we assume that the excitation operator r„ creates 
states orthogonal to |APlroG), {v\ APlroG) = 0. The 
projection manifold {zz} will depend on the choice of the 
cluster operator T (vide infra). 

The energy can be calculated by projecting against 
the reference determinant of |APlroG), i.e., multiplying 
eq. 0 by (4>o| and using intermediate normalization, 

<$ 0 |(A +[A, T])| APlroG} = E. (9) 

The only constraint on the cluster operator we have 
made so far is that it creates states that are orthogonal 
to the APlroG reference function. A possible choice for 
the cluster operator that ensures this orthogonality con¬ 
dition is to include substitutions between the occupied 
and virtual orbitals with respect to |APlroG). If only 
double excitations are included, the cluster operator is 
specified as 

^ occ virt 

^ = ( 10 ) 

ij ab 

where E a i = a) a ai + a|cq is the singlet excitation oper¬ 
ator and the cluster amplitudes are symmetric with re¬ 
spect to pair-exchange, i.e., f“ b = f b “. The prime in the 
above summations indicates that pair-excited determi¬ 
nants are excluded in the cluster operator, i.e., = 0 

(as those excitations do not fulfill the orthogonality con¬ 
dition, (3 APlroG) = 0). 
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To arrive at a computationally feasible model, we will 
further restrict the cluster operator of eq. (101 to allow 
for excitations with respect to the reference determinant 
only. Thereby, we exclude possible redundancies in ex¬ 
citations and amplitudes. The projection manifold then 
contains all doubly-excited determinants with respect to 
|$o)- Furthermore, as basis for the bra states of the pro¬ 
jection manifold, we will use the convenient choice 


lab I _ ^/abl 4. A/ ab | 
\*i I 3 I " r g 'I > 


( 11 ) 


where (“j’l = ($o| Ej b Ei a . The bra basis then forms a 
biorthogonal basis which satisfy the normalization con¬ 
dition 


iijlkl} — 3iajb,kcld "P 3jbia,kcld- (12) 

The doubles amplitudes are obtained by solving 

a linear set of equations 


B/j. + A^ v t v — 0, (13) 

V 


where the sum runs over all double excitations (with¬ 
out pair excitations) and = Bj bkc = (p-lH |APlroG), 
while A^ v = Aj bkCtiaof = l(^|[F,£' ai ^/ 0 ]|APlroG). 

The energy correction Ecorl with respect to the APlroG 
reference wavefunction is given as 

E corl = (i k \\ bc ) + (jk\bc)). (14) 

jbkc 

where we have used the standard notation for the ex¬ 
change intergrals, (jk\\bc) = (jk\bc ) — (jk\cb), and physi¬ 
cists’ notation for the two-electron integrals. 

Similarly, the contribution of single excitations can be 
accounted for by including 

f 1= Y,t“E ai (15) 

ia 

in the cluster operator. Restricting the single excitations 
to the APlroG reference determinant |$o), the singles 
projection manifold contains all singly-excited determi¬ 
nants with respect to |$o). In analogy to the doubles 
projection manifold, the bra states of the singles projec¬ 
tion manifold are chosen to form a biorthogonal basis 
with the convenient normalization condition 


(i lj) — $ai,bj I 


(16) 


where (“| = |(“| = |($ 0 | E ia . The single and double am¬ 
plitudes are obtained by solving a coupled set of linear 
equations equivalent to eq. (13 ) where fi and v now run 


over all single and double excitations. The energy cor¬ 
rection with respect to the APlroG reference value is as 
follows 

E& - 2 £ Fj b t b j + J2 t%((Mbc) + (.jk\bc )), (17) 

jb jbkc 


where Fj b are the elements of the Fock matrix, Fj b = 
hjb + Em C ((MIJ TO ) + (bm\jm)) and hj b are the one- 
electron integrals. Note that, in contrast to canonical 
Hartree-Fock orbitals, the Fock matrix is not diagonal 
when the orbitals are optimized within APlroG. In the 
APlroG-LCC approach, the single excitations thus con¬ 
tribute both directly to the energy correction and indi¬ 
rectly through coupling to the doubles equations. 

We will abbreviate the LCC correction using T = T 2 
as APlroG-LCCD, while APlroG-LCCSD indicates that 
the cluster operator contains single and double excita¬ 
tions, T = Tj -f T 2 . 


III. COMPUTATIONAL DETAILS 

A. APlroG 

All geminal calculations have been performed in the 
Horton 2.0.0 software package [56]. All restricted (vari- 
ationally) orbital-optimized APlroG calculations were 
allowed to freely relax without any spatial symmetry con¬ 
straints, Le., no point group symmetry was imposed. For 
all molecules, all orbitals were active. In the following, we 
will abbreviate (variationally) orbital-optimized APlroG 
as APlroG. 


B. PTa and PTb 

The PTa and PTb calculations were performed us¬ 
ing the Horton 2.0.0 software package [56]. For all 
molecules, the optimized APlroG natural orbitals were 
taken as orbital basis and all electrons and orbitals have 
been included in the active space. 


C. LCCD and LCCSD 


The Linearized Coupled Cluster models with double 
and single and double excitations have been implemented 
in a developer version of Horton 2.0.0 [56]. The nat¬ 
ural orbitals of APlroG were chosen as orbital basis for 
all molecules studied. Furthermore, all electrons and or¬ 
bitals were correlated. 


D. Coupled Cluster 

The Coupled Cluster Doubles (CCD), CC Singles and 
Doubles (CCSD) as well as CC Singles, Doubles and per¬ 
turbative Triples (CCSD(T)) calculations have been car¬ 
ried out in the DALTON2013 software package (57]. In 
each case, all electrons and orbitals were correlated and 
no spatial symmetry was imposed. 
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E. Relativity and Basis Sets 

For the C 2 and F 2 molecules, Dunning’s aug-cc-pVDZ 
(C, F:(10s5p2d) —> [4s3p2d]) and aug-cc-pVTZ (C, F: 
(Ils6p3d2f) —► [5s4p3d2f]) basis sets were used, while 
the STO- 6 G basis set was used for the H atoms in H 50 
to allow for a comparison to DMRG reference data. 

For UOj" 1 ", scalar relativistic effects were incorporated 
through relativistic effective core potentials (RECP). In 
all calculations, we have used a small core (SC) RECP 
(60 electrons in the core) with the following contraction 
scheme (12sllpl0d8f) —► [8s7p6d4f]. For the lighter el¬ 
ements (O), the cc-pVDZ basis set of Dunning was em¬ 
ployed, (10s5pld) -A [4s3pld]. 

F. Fitting procedure 

The potential energy curves of diatomic molecules were 
obtained by varying bond lengths in a range of 1 . 2 —4.0 A 
and 1.1 —3.2 A for the Ft and C 2 molecules, respectively. 
The points on the resulting potential energy curve were 
used for a subsequent generalized Morse function [ 58] fit 
to obtain the equilibrium bond lengths (R e ) and potential 
energy depths (D e ). The harmonic vibrational frequen¬ 
cies (w e ) were calculated numerically using the five-point 
finite difference stencil [59]. 

IV. NUMERICAL RESULTS 
A. Dissociation of C2 

The carbon dimer is one of the most complex diatomic 
molecules that can be formed from the first-row elements 
of the periodic table. The unusual nature of the carbon- 
carbon bond and the question concerning its bond order 
attracted a lot of attention from theoretical chemists EH- 
[541 in recent years. Around the equilibrium distance, the 



Distance [A] 

FIG. 1. Potential energy surfaces for the dissociation of the 
C 2 molecule using an aug-cc-pVTZ basis set. 


electronic structure of the C 2 molecule has two dominant 
configurations lagla^2ag2a^lir^ and la^lcrl2agln^Sa'g 
as well as a number of other configurations arising from 
low-lying excited states ESI E5HS3- When the two 
carbon atoms are pulled apart, the molecular system 
becomes strongly multi-reference. However, even for 
stretched carbon-carbon distances, dynamical electron 
correlation effects remain non-negligible B51E51EZI- a 
reliable theoretical description of spectroscopic constants 
(bond lengths, potential energy well depths, and vibra¬ 
tional frequencies) thus requires a balanced treatment of 
all types of electron correlation effects |441155] (static, 
non-dynamic, and dynamic) along the dissociation path¬ 
way. Since highly accurate reference data for the spectro¬ 
scopic constants of C 2 is available, it is an ideal candidate 
to test our APlroG-LCC approach. 

The upper part of Table [I] summarizes the spectro¬ 
scopic constants of the C 2 molecule for different basis 
sets and quantum chemistry methods including various 
dynamical correlation models based on an APlroG refer¬ 
ence function. As expected, APlroG predicts too short 
equilibrium bond distances and too shallow potential well 
depths for all basis sets studied. If dynamical correla¬ 
tion is included a posteriori on top of the APlroG ref¬ 
erence function using perturbation theory, spectroscopic 
constants improve only slightly compared to MRCI-SD 
reference data. Although APlroG-PTa predicts equilib¬ 
rium bond distances that are in perfect agreement with 
MRCI-SD, the potential well depth is overestimated and 
the differences with respect to MRCI-SD are even larger 
than for APlroG without PTa correction. Furthermore, 
APlroG-PTb does not improve potential well depths and 
vibrational frequencies compared to APlroG when the 
basis set is enlarged from aug-cc-pVDZ to aug-cc-pVTZ. 
In contrast to the PTa and PTb models, an LCC correc¬ 
tion on top of APlroG including doubles and singles and 
doubles yields spectroscopic constants that are in very 
good agreement with MRCI-SD reference data, outper¬ 
forming NEVPT2 (differences are less than 2 kcal/mol 
for potential well depths using aug-cc-pVTZ). 

Figure [l] shows the fitted potential energy surfaces for 
the aug-cc-pVTZ basis set. We should note that all po¬ 
tential energy surfaces were adjusted to Zero in the dis¬ 
sociation limit. Comparing the spectroscopic constants 
of Table |T] we can conclude that the MRCI-SD poten¬ 
tial energy surface would lie between the APlroG-LCCD 
and APlroG-LCCSD potential energy curves. All other 
quantum chemistry methods yield potential energy sur¬ 
faces that deviate more from the expected MRCI-SD ref¬ 
erence curve. 


B. Dissociation of F2 

F 2 is a well-known example of diatomic molecules 
where dynamical electron correlation effects play a dom¬ 
inant role [M]. Furthermore, theoretical studies indicate 
that large basis sets and robust dynamical electron cor- 





















5 


TABLE I. Spectroscopic constants for the dissociation of the C 2 and F 2 molecule for different quantum chemistry methods and 
basis sets. The differences are with respect to MRCI-SD reference data EDI- 



Method 

aug-cc-pVDZ 

MAl d* MI 

oo e [cm x ] 

aug-cc-pVTZ 

MAI[^f] 

00 e [cm 

c 2 

APlroG 

APlroG-PTa 

APlroG-PTb 

APlroG-LCCD 

APlroG-LCCSD 

NEVPT2 

1.240(—0.033) 104.6(—15.7) 1943(+136) 
1.273(+0.000) 152.7(+22.4) 1902(+95) 
1.260(—0.013) 116.3(—14.0) 1863(+56) 
1.261(—0.012) 124.3(—7.0) 1926(+119) 
1.266(—0.007) 134.5(+4.2) 1855(+48) 
1.259(—0.014) 135.0(+4.7) 1924(+117) 

1.227(—0.025) 132.9(—8.2) 1780(-56) 
1.251(+0.001) 160.4(+19.3) 1889(+53) 
1.235(—0.017) 127.2(—13.9) 1938(+102) 
1.240(—0.012) 139.3(—1.8) 1916(+80) 
1.240(—0.012) 143.0(+1.9) 1926(+90) 
1.244(—0.008) 148.0(+6.9) 1886(+50) 


MRCI-SD 

1.273 

130.3 

1807 

1.252 

141.1 

1836 

f 2 

APlroG 

APlroG-PTa 

APlroG-PTb 

APlroG-LCCD 

APlroG-LCCSD 

CCSD 

CCSD(T) 

1.521(+0.068) 
1.398(—0.055) 
1.473(+0.020) 
1.466(+0.013) 
1.462(+0.009) 
1.426(—0.017) 
1.450(—0.003) 

12.8(—15.7) 
30.1(+1.6) 

39.5(+11.0) 

40.1(+11.6) 

57.3(+18.8) 

886(+85) 
636(—165) 
832(+31) 
780(—21) 
793(—8) 
917(+116) 
841(+40) 

1.467(+0.047) 
1.448(+0.028) 
1.417(—0.003) 
1.433(+0.013) 
1.431(+0.011) 
1.396(—0.024) 
1.419(—0.001) 

16.2(—7.7) 
33.7(—0.2) 

45.5(+11.6) 

46.7(+12.8) 

69.2(+35.3) 

703(—189) 
847(—45) 
891(—1) 
872(—20) 
883(—9) 
1004(+12) 
911(+19) 


MRCI-SD 

1.453 

28.5 

801 

1.420 

33.9 

892 


relation models are required to reproduce the experimen¬ 
tally determined spectroscopic constants |69lf7T| . These 
features make the F 2 molecule a good test case to assess 
the reliability of the LCC correction on top of an APlroG 
reference function. 

The bottom part of Table |T] summarizes the spectro¬ 
scopic constants for the dissociation process of the F 2 
molecule using different basis sets and dynamical corre¬ 
lation models. We should note that both CCSD(T) and 
APlroG-PTb diverge in the dissociation limit and thus 
the potential energy well depths are not given in Table |T] 
(see also Figure [2j|. To obtain an estimate for D e , we 
have taken the MRCI-SD results by Peterson |60| . Note 
that r e and oj e as predicted by MRCI-SD are in good 
agreement with CCSD(T) calculations. 

As observed for the C 2 molecule, APlroG consider¬ 
ably overestimates the equilibrium bond length and un¬ 
derestimates the potential energy depth, which can be 
attributed to the large fraction of dynamical correlation 
that cannot be captured by restricting the wavefunction 
to electron-pair states. Although APlroG-PTa yields po¬ 
tential energy depths that are in very good agreement 
with MRCI-SD reference data, equilibrium bond lengths 
and vibrational frequencies deviate considerably from the 
MRCI-SD reference values. Specifically, when increas¬ 
ing the basis set from aug-cc-pVDZ to aug-cc-pVTZ, the 
predicted equilibrium bond length changes from being 
too short to being overestimated. In contrast to PTa, 
APlroG-PTb results in equilibrium bond lengths and vi¬ 
brational frequencies that are in perfect agreement with 
MRCI-SD reference data when the basis set is increased 
to triple-zeta quality, fails however in the vicinity of dis¬ 
sociation. The LCC correction on top of APlroG re¬ 
sults in the most stable and robust dynamical correla¬ 
tion model, yielding similar results for increasing basis 
set sizes and outperforming CCSD. Specifically, equilib- 



Distance [A] 

FIG. 2. Potential energy surfaces for the dissociation of the 
F 2 molecule using a aug-cc-pVTZ basis set. 


rium bond lengths and vibrational frequencies are in very 
good agreement with MRCI-SD reference data. However, 
potential energy well depths are considerably overesti¬ 
mated by APlroG-LCC (around 12 kcal/mol with re¬ 
spect to MRCI-SD), which might be attributed to re¬ 
stricting the cluster operator to singles and doubles ex¬ 
citations [69| (cf ., CCSD overestimates D e by more than 
20-35 kcal/mol, depending on the basis set size). 


C. Symmetric Dissociation of H50 

The symmetric stretching of hydrogen chains is com¬ 
monly used as a molecular model for strongly correlated 
systems and remains a challenging problem for conven¬ 
tional quantum-chemistry methods [5711321150-93 . Re¬ 
cently, we have shown that APlroG accurately describes 
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FIG. 3. Symmetric dissociation of the H 50 chain using the 
STO- 6 G basis set obtained from different methods. The 
DMRG reference data are taken from Ref. [90) . while the MP2 
and APlroG data are taken from Ref. EZl- 
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the potential energy surface of the symmetric dissocia¬ 
tion of H 50 in the vicinity of dissociation, but deviates 
from DMRG reference data around the equilibrium and 
for stretched interatomic distances mi which can be at¬ 
tributed to the missing dynamical correlation energy. 

Table [Tl] summarizes the non-parallelity error (NPE) 
per hydrogen atom for the symmetric dissociation of H 50 
obtained by APlroG and different dynamical correlation 
models. The large NPE per hydrogen atom of APlroG 
can be associated with the missing dynamical correlation 
energy around the equilibrium distance and for stretched 
interatomic bond lengths. Adding dynamical correlation 
a posteriori improves the NPE per hydrogen atom con¬ 
siderably. While PTa, PTb, and LCCD have a similar 
NPE per hydrogen atom of about 4.5 m Eh, the NPE 
is reduced to less than 1.5 m Eh if single excitations are 
included in the cluster operator. 

The importance of single excitations in the LCC model 
is also noticeable in the shape of the potential energy sur¬ 
face shown in Figure [3j While APlroG-PTa, APlroG- 
PTb, and APlroG-LCCD have similar potential energy 
curves (in terms of shape and total electronic energies) 
with APlroG-PTb being lowest in energy, the potential 
energy surface obtained by APlroG-LCCSD considerably 
deviates from the aforementioned dynamical correlation 
models for short and intermediate bond lengths (around 
1.0 and 1.5 A). Furthermore, the APlroG-LCCSD poten¬ 
tial energy curve is in very good agreement with DMRG 
reference data (see Figure [3]). 


D. Symmetric dissociation of UO^ 

The uranyl cation (UO^) is a small building block of 
a large variety of uranium-containing complexes EMzzj. 
This molecule has a linear structure and a singlet ground- 
state electronic configuration. Its characteristic symmet- 


TABLE II. Non-parallelity error per hydrogen atom for the 
symmetric dissociation of the H 50 chain with respect to 
DMRG reference data. The MP2 and APlroG data are taken 
from Ref. [37]. 


Method 

NPE/H [m E h ] 

MP2 

43.540 

APlroG 

6.187 

APlroG-PTa 

4.831 

APlroG-PTb 

4.170 

APlroG-LCCD 

4.506 

APlroG-LCCSD 

1.389 


TABLE III. Spectroscopic constants for the symmetric disso¬ 
ciation of the UOj” 1 " molecule for different quantum chemistry. 
The differences are with respect to CCSD(T) reference data. 
The CASSCF and CC data were taken from Ref. 1461 


Method 


w 


APlroG 

APlroG-PTb 

APlroG-LCCD 

CAS(10,10)SCF 

CAS(12,12)SCF 

CCD 

CCSD 


1.669(—0.047) 
1.715(—0.001) 
1.708(—0.008) 
1.694(—0.022) 
1.707(—0.009) 
1.690(—0.026) 
1.697(—0.019) 


1062(+53) 

1340(+331) 

1073(+64) 

1079(+70) 

1034(+25) 

1125(+116) 

1068(+59) 


CCSD(T) 


1.716 


1009 


ric and asymmetric U-0 vibrational frequencies are used 
to identify the presence of UO^ in larger molecular as¬ 
semblies IZSEHIESI- While the electronic structure of the 
uranyl cation is well-known around the equilibrium struc¬ 
ture [36] I50H85] , the complicated nature of the U-0 bond 
hampers a theoretical description at larger U-0 distances 
using standard quantum chemistry approaches [46] . One 
of the limiting factors that impede theoretical studies 
is the large number of strongly-correlated electrons dis¬ 
tributed among 5/-, 6 d-, and 7s-orbitals. In addition, 
the 6 s- and 6 p-core-valence orbitals are easily polarizable 
and have a non-negligible contribution to the correlation 
energy. However, around the equilibrium structure, the 
uranyl cation is well described by single-reference CC the¬ 
ory if all important electrons are correlated. This allows 
us to assess the performance of the LCC models in de¬ 
scribing dynamical correlation effects ordinating from the 
5/-, 6 d-, and 7s- as well as the core-valence electrons. 

The equilibrium bond lengths and vibrational frequen¬ 
cies obtained by different quantum chemistry methods 
are shown in Table |III| As expected, APlroG consider¬ 
ably underestimates the equilibrium bond length, while 
uj e is in good agreement with CCSD(T). Adding dynam¬ 
ical correlation effects on top of APlroG shifts r e closer 
to CCSD(T) reference data. The shape of the poten¬ 
tial, however, strongly depends on the APlroG dynam¬ 
ical correlation model. Specifically, PTb results in a 
much steeper potential energy surface overestimating vi¬ 
brational frequencies by more than 330 cm -1 compared 
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Distance [A] 

FIG. 4. Potential energy surfaces for the symmetric stretching 
of the U02+ molecule around the equilibrium geometry. Note 
that the CASSCF potential energy surfaces are much higher 
in energy and are thus not shown. 


to CCSD(T), while LCCD preserves the shape of the po¬ 
tential energy surface and yields a vibrational frequency 
that agrees well with APlroG and CCSD(T) data (dif¬ 
ferences amount to ca. 60 cm -1 ). The overall accuracy 
of APlroG-LCCD lies between CCSD and CCSD(T), be¬ 
ing closer to the latter. We should emphasize that PTa 
completely fails for the UO 2 molecule and produces a 
discontinuous potential energy surface around the equi¬ 
librium (see also Figure [4]). Furthermore, the CASSCF 
equilibrium distance strongly depends on the size of the 
active space chosen in CASSCF calculations. Specifi¬ 
cally, increasing the active space from CAS(10,10) to 
CAS(12,12), i.e., including the a- and ( 7 *-orbitals, results 
in spectroscopic constants that are in good agreement 
with APlroG-LCCD and CCSD(T) data. 

Figure [4] shows the fitted potential energy surfaces 
around the equilibrium for selected quantum chemistry 
methods. APlroG-LCCD yields total electronic energies 
that are between CCSD and CCSD(T), while the poten¬ 
tial energy surface predicted by AP IroG-PTb is consider¬ 
ably lower than the CCSD(T) reference curve. Note that 
the potential energy surfaces optimized by CASSCF lie 
much higher in energy and are thus not shown in Fig¬ 
ure m 


V. CONCLUSIONS 

We have presented an alternative model to capture dy¬ 
namical correlation effects on top of an APlroG reference 
functions that uses a Linearized Coupled Cluster ansatz. 
Our approach is motivated by previous studies of an LCC 
correction to an APSG reference function [S3]. Specifi¬ 
cally, our cluster operator is restricted to doubles and 
singles and doubles excitations as in the standard Cou¬ 
pled Cluster approach, i.e., excitations from the occupied 
to the virtual orbitals of some reference determinant. We 
have compared this new dynamical correlation ansatz to 
the PTa and PTb perturbation theory models as well as 
standard quantum chemistry approaches for the C 2 and 
F 2 molecules, the H 50 hydrogen chain, and the uranyl 
cation, UO^. 

In general, both APlroG-LCC models yield similar 
spectroscopic constants and are closest to MRCI-SD, 
CCSD(T), and DMRG reference data for all molecules 
we have investigated. Furthermore, LCCD and LCCSD 
represent more robust and reliable dynamical correlation 
models than PTa and PTb. Our study demonstrates that 
the success and failures of PTa and PTb are difficult to 
anticipate a priori and are strongly system-dependent. 
While PTa yields reasonable results for equilibrium bond 
lengths of the C 2 molecule and for the potential energy 
depth of the F 2 molecule, the corresponding C 2 potential 
energy depth and the F 2 equilibrium bond length signif¬ 
icantly differ from reference data. A similar behavior 
was observed for PTb. In contrast to the perturbation 
theory models, the LCC ansatz is able to capture differ¬ 
ent flavours of dynamical correlation effects reliably, as 
present in the C 2 , F 2 , H 50 , and U 02 + molecules. 
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